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Abstract. We study static cosmic string solutions in the context of Bekenstein-type 
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in Bekenstein-type models strongly depart from the standard Nielsen-Olesen solution 
^\J ' with the electromagnetic energy concentrated along the string core seeding spatial 

variations of the fine structure constant, a. We consider models with a generic gauge 
vj I kinetic function and show that equivalence principle constraints impose tight limits on 

^vq ' the allowed variations of a induced by string networks on cosmological scales. 
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1. Introduction 

Cosmological theories motivated by models with extra-spatial dimensions ^Ij have 
recently attracted much attention, in large measure due to the possibility of space- 
time variations of the so-called "constants" of Nature J2l 13 IH IH]- Among these are 
varying-a models (where a is the fine-structure constant) such as the one proposed by 
Bekenstein [H]. In Bekenstein-type models the variation of the fine structure constant, 
a = e^/47r h c, is sourced by a scalar field, if, coupling minimally to the metric and to the 
electromagnetic term, "_p2" ^ ^j ^ non-trivial gauge kinetic function, Bp{ip). The interest 
in this type of models has been increased with recent results coming from both quasar 
absorption systems [71 |H] (see however P 1101) and the Oklo natural nuclear reactor 
[11 suggesting a cosmological variation of a at low red-shifts. Other constraints at low 
redshift include atomic clocks ^2] and meteorites ^Hj- At high redshifts there are also 
upper limits to the allowed variations of a coming from either the Cosmic Microwave 
Background or Big Bang Nucleossynthesis [H US CEl El CHI CHI EH • 

Meanwhile on the theoretical side some effort has been made on the construction 
of models which can explain a variation of a at redshifts z ~ 1 — 3 of the same order as 
that reported in refs. [3 JH] while being consistent with the other constraints at lower 
redshifts. Although this cannot be achieved in the simplest class of models in which 
the potential and gauge kinetic function are expanded up to linear order |^ there are 
other classes of models which seem in better agreement with the data (see for example 
122] )■ Moreover, in some of these models |2Ill221l2nil21l2Sll2Sll2ZI the fine structure 
constant is directly related to the scalar field responsible for the recent acceleration of 
the Universe ^11^ . 

A number of authors have also studied the spatial variations of the fine structure 
constant induced by fiuctuations in the matter fields showing that they are proportional 
to the gravitational potential and are typically very small to be detected directly with 
present day technology except perhaps in the vicinity of compact objects with strong 
gravitational fields (see for example jHOl l3T | l32 l l33] ) . 

The variation of a in the early Universe has even been associated with the solutions 
to some of the problems of the standard cosmological model jSU JHSl- However, it has 
been shown that the ability of specific models to solve some of the problems of the 
standard cosmology (in particular the horizon and fiatness problems) is directly related 
to the evolution of 'cosmic numbers' which are dimensionless parameters involving 
cosmological quantities rather than the evolution of dimensionless combinations of the 
so-called 'fundamental constants of nature' [86ll87j. 

Note that a change in the value of the fine structure constant should in principle 
be accompanied by a variation of other fundamental constants as well as the grand 
unification scale. However, we shall adopt a phenomenological approach and neglect 
possible variations of other fundamental constants. 

In this article we consider the case of static Abelian vortex solutions whose 
electromagnetic energy is localized along a stable string-like core which acts as a source 
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for spatial variations of a in the vicinity of the string. We generahze to a generic 
gauge kinetic function the work of ref. |HH1 and study the hmits imposed by the Weak 
Equivalence Principle p^ W)\ on the allowed cosmological variations of a. The article 
is organized as follows. In Sec. II we briefly introduce Bekenstein-type models and 
obtain the equations describing a static string solution. We describe the numerical 
results obtained for a number of possible choices of the gauge kinetic function in Sec. 
Ill discussing the possible cosmological implications of string networks of this type in 
the light of equivalence principle constraints in Sec. IV. Finally, in Sec. V we summarize 
our results and briefly discuss further prospects. Throughout this paper we shall use 
units in which h = c = 1 and the metric signature H . 

2. Bekenstein-type models 

We first review Bekenstein-type models with a charged complex scalar field (j) undergoing 
spontaneous symmetry breaking. Let be a complex scalar field with a U(l) gauge 
symmetry and a^ be the gauge field. Let us also assume that the electric charge is a 
function of space and time coordinates, e = eoe{x^) where e is a real scalar field and cq 
is an arbitrary constant charge. The Lagrangian density in Bekenstein-type models is 
given by (see for example |SH|): 

c = {D,<pnD'^<p)-v{<p) -^^Bp{v)f,,r 
+ \d,^d^^, (1) 

where Bp^Lp) = e^'^i^p) is the gauge kinetic function and yj is a scalar field. In eqn. (|H), 
D^ (j) = {d^ — ieoo^) (f) are covariant derivatives and the electromagnetic field tensor is 
given by 

f^lu = df,a^ - duO^ . (2) 

The function Bp = e~^(v^) acts as the effective dielectric permittivity which can be 
phenomenologically taken to be an arbitrary function of if. We assume that V{(f)) is the 
usual Mexican hat potential with 

V{<P) = \(\^\' + ^] , (3) 



where A > and m? < are constant parameters and i] = mj \/2\ is the symmetry 
breaking scale. The Lagrangian density in eqn. (^ is then invariant under U(l) gauge 
transformations of the form 5 0= —iC, 0, Op, ^> a^.^ C,^p,. 

Zero variation of the action with respect to the complex conjugate of 0, i.e., 0*, 
gives: 

B.B'^=-|^- (4) 

Variation with respect to a^ leads to: 

d^Bp{^)r\ = J^ (5) 
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with the current j^ defined as 

Finally, variation with respect to (f gives 



1 dBFJip ) ^2 



(6) 



(7) 



We now look for the static vortex solutions in these theories adopting the following 
ansatz 

<P = X{r) e^"' , (8) 

ag = a{r) , (9) 

with all other components of a^ set to zero. Here we are using cylindrical coordinates 
{r,9,z), n is the winding number and x(r) is a real function of r. Substituting the 
ansatz given in (jHHHl) into eqns. (jHEj) one gets 



1 d 
r dr 



,dx 
dr ^ 

1 d 



Co a] + rri^ + 2Ax^ 



X = 0, 



^(BF-—{ra) 
dr \ r dr 



n 



+ 2eo ( - - eoa ) X 



I d ( d'^\ IdBpi^p) fl d 
r dr \ dr J 2 dip \r dr 



0. 



Note that 



r = 2r /, 



rd 



1 d 
r dr 



1 2 



ra 



(10) 

(11) 
(12) 

(13) 



which substituted in eqn. ((Tj) gives the factor of one-half in the second term on the left 
hand side of eqn. ()12|1 . This corrects an error in ref. |2Si which neglected the factor of 
2 in eqn. (fTT^ . 

We also investigate the dependence of the energy density on the radial coordinate r. 
For static strings the stress-energy tensor take a diagonal form with the energy density 
of the vortex being p = Tq = q^'^Tqq = Tqq, with 

/,\2 /,\2 , /,\2 

ld\) 
P 



dr I \dr J 



2 r^ \dr , 



+ 



n 



eQV\ 



■) 



X' + A X' + 



m 
2A 



(14) 



while the spatial components of the stress-energy tensor are given by T- = diag{pr, pe, Pz) 
withp^ = —p (here v = ar). Therefore, the energy density of the vortex, p, is everywhere 
positive, while the longitudinal pressure pz is negative. In fact, this is also one of the 
defining features of canonical cosmic strings. 
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3. Numerical Solutions 

We consider solutions to the coupled non-linear equations for a static straight string. 
Since no exact analytic solution has yet been found it proves useful to reduce equations 
()10II12|) to a set of first order differential equations for numerical implementation. 
Let us introduce three new variables 

^=&r, (16) 

ar 

ar 
with V = ar. Then by substituting p5lll7j) into equations ()1()II12J] . one gets 

X, (18) 



da o 

dr r 

db _ 1 
~dr ~ ^ 



' + m +2 Ax 



r]b-2eo[ x 



dip 



(19) 



di] 7] 1 dBp 2 

dr r 2 dip 

Then we have, at all, a set of six ordinary first order differential equations, which 
requires, at least, six boundary conditions to be solved numerically. The appropriate 
boundary conditions are [IH] 

(21) 

(22) 

(23) 

Therefore we have a two point boundary value problem with four conditions at the 
origin and two conditions far from the core. Let us start by discussing the boundary 
conditions far from the core (r -^ +oo). In this limit the scalar field has a constant 



limx(r) = 




}^^^\Kr) -\^ ^^ 


limw(r) = 


lim v[r) = — , 
eo 


limSir(r) = 1 


lim ri = . 



value given by x = v ~'^^/(2A) and from equation (|T(Hl one immediately sees that this 
implies that v{r) = n/eo far from the core. There are also four boundary conditions at 
the string core. In this limit, the phase 6 in equation (jH)) is undefined which implies that 
X must vanish at r = 0. Also, v must vanish at the string core. Otherwise the magnetic 
energy density would diverge at the core of the string. We also normalize the electric 
charge such that e = Cq at the string core (this gives the boundary condition for Bp). 
Finally, the boundary condition for rj becomes evident by using the Gauss law to solve 
equation (fT^ assuming that there are no sources of a variation other than the string. 

In order to solve this problem numerically we used the relaxation method which 
replaces the set of six ordinary differential equations by finite-difference equations on a 
mesh of points covering the range of the integration. This method is very efficient if a 
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Figure 1. The numerical solution of the scalar field ?/> = In e as a function of distance, 
r, to the core of string, for the original Bekenstein model. If w < then e — > cx) when 
r — > oo. On the other hand, if a; > then e ^ when r — > oo. The dashed line 
represents the constant-a theory, which corresponds to the limit uj — > c». 



good initial guess is supplied. In our case, the solutions of the standard Nielsen-Olesen 
vortex can be used to generate a good initial guess. We checked that our code reproduces 
the results for the standard Nielsen-Olesen vortex if Bp = 1. Throughout this paper we 
shall take A = 1/2 and n = 1 for definiteness and use units in which h = c = —m? = 1. 

3.1. Exponential coupling 

The general prescription detailed above can be particularized to specific choices of gauge 
kinetic functions. First, let us consider 

Bpi^) = e-^ . (24) 

Then the Lagrangian density (0) becomes 

1 ff^" f 
C= (D^^r {D,<P)-V{<P) --L-p^ + 

+ ^d,ed^e, (25) 

which is now written in terms of the field 

e = - = exp (^] (26) 

and recovers the original Bekenstein model (studied in detail by Magueijo et al in ref. 
pS] in a similar context). In this model a; is a coupling constant. 

It is easy to show that in the limit u; — > oo one recovers the Nielsen-Olesen vortex 
with constant e. Although the gauge kinetic function in eqn. (j^^ is only well defined 
for a; > the model described by the Lagrangian density in eqn. (J25|l allows for both 
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negative and positive values of u. However, note that if u; < the energy density is 
no longer positive definite. In Fig^ we plot the numerical solution of the scalar field 
ip = lne as a. function of distance, r, to the core of string, in the context of the original 
Bekenstein model. Note that ii u < then e diverges asymptotically away from the 
string core. On the other hand if a; > then ip goes to zero when r — * oo. In the large 
u limit the curves for positive and negative u are nearly symmetric approaching the 
dashed line representing the constant-a model when a; ^ oo. 

3.2. Polynomial coupling 

Another example of a class of gauge kinetic functions is given by 

N 

Bp{^) = 1.0 + J2P,^\ (27) 

j=i 

in which j3i are dimensionless coupling constants and A^ is an integer. If /5i = it 

is easy to verify that the classical Nielsen- Olesen vortex solution with constant a is 

still a valid solution. This means that there is a class of gauge kinetic functions for 

which the classical static solution is maintained despite the modifications to the model. 

Substituting the gauge kinetic function in eqn. (fT^ one gets 



1 d dip 
r dr dr 




N \ 

k=l J 

N \ 

^(2A;)/32.^2.-i _ (28) 

fc=i / 

which shows that the transformation /?j -^ —f3i for odd i modifies the sign of the solution 
of ip{r) without changing x or 6 since Bp is kept invariant. 

We will see that both for /3i > and /5i < the behaviour of ■?/' = In e is similar 
to that of the original Bekenstein model described by the lagrangian density in eqn. 
fl25|l . with u; > 0, in particular in the limit of small |/?i| /large uo. In fact a polynomial 
expansion of the exponential gauge kinetic function of the original Bekenstein model 
(see eqn. (I24|l ) has (3i = —2/-s/lj. This relation between the models arises from the 
fact that the exponential coupling by Bekenstein theory can be expanded in a series of 
powers of ip according to 

-2y 

On the other hand, as mentioned before, if /3i = one recovers the standard result 
for Nielsen-Olesen vortex with constant-a, independently of the chosen values of f3i for 
i > 1. 

We have studied the behaviour of the solutions of eqns. (|15II2(J|) for various values 
of A^ but for simplicity we shall only consider A^ < 2 in this paper. In particular, we 
consider 

Bp = l.O + (3i^ + (32¥>^ (30) 



A = ^ • (29) 
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Figure 2. The numerical solution of the scalar field ^ = In e as a function of distance, 
r, to the core of string, for a polynomial gauge kinetic function. Models 0, 1 and 
2 are defined by /3i — {(32 arbitrary), /3i = — 3, /32 = (linear coupling) and 
Pi = —5, /32 = 10 respectively. 




Figure 3. The numerical solution of the fields x{^) a-^d b{r) as a function of distance, 
r, to the core of string, for models 0, 1 and 2. Note that the change in b{r) with respect 
to the standard constant-a result is much more dramatic than the change in %('')• 
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Figure 4. Plot of Cii''^) = log(Xi/Xo) for the different polynomial gauge kinetic 
functions. Note that although this was not very visible in Fig. |3 even a small value 
of /3i leads to a change of the vortex solution. 



with two free parameters. 

In Fig. 121 we plot the numerical solution of the scalar field ip = \ne as a function 
of distance, r, to the core of string, for a polynomial gauge kinetic function. Models 1 
and 2 are defined Pi = —3, /92 = (linear coupling) and Pi = —5, P2 = 10 respectively. 
Model (dashed line) represents any model with (3i = and has a = constant. Note 
that the replacement (3i -^ —f3i does not modify the solution for ip. 

In Fig. IHlwe plot the numerical solution of the fields x(^) and b{r) as a function of 
distance, r, to the core of string, for models 0, 1 and 2. We see that the change in 6(r) 
with respect to the standard constant-a result is much more dramatic than the change 
in x(^)- In order to verify the modification to x{^) in more detail we define the function 



0(r) 




(31) 



and plot in Fig. |3]the results for the different models. We clearly see that even a small 
value of Pi leads to a modification of the vortex solution with respect to the standard 
Nielsen-Olesen solution. 

Finally, we have also studied the behaviour of the energy density in this model. 
In fact, since the fine structure constant varies we have a new contribution due to the 
field ip, to the total energy of the topological defect. As has been previously discussed 
in ref. ^j the contributions to the energy density of the string can be divided into 
two components. One component is localized around the string core (the local string 
component) the other is related to the contribution of the kinetic term associated with 
the spatial variations of the fine structure constant and is not localized in the core of 
the string. The energy profile of this last contribution is analogous to that of a global 
string, whose energy per unit length diverges asymptotically far from the core. In Fig. 
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Figure 5. The energy density as a function of the distance, r to the string core for 
models 0, 1 and 2. The dashed hne represents the constant-a model These results 
clearly show an increase of the energy density due to the contribution of the extra field 

if. 

|S]we plot the string energy density as a function of the distance, r to the string core for 
models 0, 1 and 2. The dashed line represents the constant-a model. We clearly see an 
increase of the energy density due to the contribution of the extra field cp. 

As an aside, let us point out that the Lagrangian density in (jT} can be generalized 
to include a potential V{ip). If V{(p) has a minimum and it is steep enough near it then 
the spatial variations on the value of the fine structure constant can be significantly 
reduced. Also, if V{ip) has more than one minimum then the Lagrangian density in (^ 
admits domain wall solutions which separate regions with different values of the fine 
structure constant. In the absence of an electromagnetic field these solutions reduce 
to the standard ones. On the other hand, if there are sources of electromagnetic field 
present (such as cosmic strings) then there may also be small fluctuations on the the 
value of the fine structure constant generated due to the electromagnetic source term in 
eqn. ((Zj). These fluctuations are expected to be more localized due to the influence of 
the potential V{ip). 

4. Cosmological implications of varying-a strings 

Having discussed in previous sections static string solutions in the context of Bekenstein- 
type models it is now interesting to investigate if such cosmic string networks can induce 
measurable space-time variations of a in a cosmological setting. In order to answer this 
question we recall that (f satisfies the Poisson equation: 

VV=^/3f. (32) 
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In eqn. (|32j) we have assumed for simplicity that the gauge kinetic function is a hnear 
function of ip with Bp{ip) = 1 + (3ip. Hence the variation of the fine structure constant 
away from the string core is given by 

27rr^ = /5/(r)/i(r^aJ , (33) 

dr 



where 



and 



^i{r) = 27r / p{r')r'dr' , (34) 







Hr] "" 



, . ^ rnrYdr', (35) 

is a function of r smaller than unity. Here Tmax represents a cut-off scale which is in a 
cosmological context of the order of the string correlation length. Far away from the 
string core /(r) is a slowly varying function of r which is always smaller than unity. 
An approximate solution for the behaviour of the field (p may be obtained by taking 
/(r) ~ const 

^^Mr^,„(M, (36) 

2n \ro/ 

where tq is an integration constant. Since e = Bf{^)~^^'^ we have 

e^l- ^'^/^(^-) In (L) , (37) 

We see that the variation of the fine structure constant away from the string core is 
proportional to the gravitational potential induced by the strings. The value of Gfi is 
constrained to be small (< 10^^ — 10^^) in order to avoid conflict with CMB and LSS 
results 1^ 1321 HSl HH m] (or even smaller depending on the decaying channels available 
to the cosmic string network |3E])- Note that these constraints are for standard local 
strings described by the Nambu-Goto action. Even though the cosmic strings studied 
in our paper are non-standard, having a local and a global component, we expect the 
limits on G/i(rmax) to be similar as those on Gfi for standard local strings. On the 
other hand the factor (3 appearing in eqn. (jHTj) is constrained by equivalence principle 
tests to be \/3\ < 10"^ G^/^ jSHl EHl ED]- Hence, taking into account that we cannot 
observe scales larger than the horizon (~ 10'* Mpc) and are unlikely to probe variations 
of a at a distance much smaller than 1 pc from a cosmic string, eqn. (jHTj) implies that 
a conservative overall limit on observable variations of a seeded by cosmic strings is 
Aa/a < 10~^^ which is too small to have any significant cosmological impact. We 
thus see that, even allowing for a large contribution coming from the logarithmic factor 
in eqn. (jTFj) . the spatial variations of a induced by such strings are too small to be 
detectable. In particular, this result means that despite the claim to the contrary in 
ref. |j2Sj cosmic strings will be unable to generate inhomogeneities in a which could be 
responsible for the generation of inhomogeneous reionization scenarios. 
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5. Conclusion 

In this article we studied Nielsen- Olesen vortex solutions in Bekenstein-type models 
considering models with a generic gauge kinetic function. We showed that there is 
a class of models of this type for which the classical Nielsen- Olesen vortex is still a 
valid solution (with no a variation). However, in general, spatial variations of a will 
be sourced by the electromagnetic energy concentrated along the string core. These 
are roughly proportional to the gravitational potential induced by the strings which is 
constrained to be small. We have shown that Equivalence Principle constraints impose 
tight limits on the allowed variations of a on cosmological scales induced by cosmic 
string networks of this type. As mentioned in ref. |SH] other defect solutions may be 
considered such as monopoles or textures but in this case we expect similar conclusions 
to be drawn as far as the cosmological consequences of the theories are concerned. 
However, despite the claim in ref. [3H] that domain walls can not be associated with 
changing-alpha theories if the scalar field ip is endowed with a potential with a Z2 
symmetry then domain-walls may form and therefore large scale inhomogeneities in the 
value of a associated with different domains may be generated. We shall leave for future 
work a more detailed study of other defect solutions and cosmological implications. 
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